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A GENERALISED GAUSS CIRCLE PROBLEM AND 
INTEGRATED DENSITY OF STATES 

JEAN LAGACE & LEONID PARNOVSKI 


Abstract. Counting lattice points inside a ball of large radius in Euclidean 
space is a classical problem in analytic number theory, dating back to Gauss. 
We propose a variation on this problem: studying the asymptotics of the 
measure of an integer lattice of affine planes inside a ball. The first term is 
the volume of the ball; we study the size of the remainder term. While the 
classical problem is equivalent to counting eigenvalues of the Laplace operator 
on the torus, our variation corresponds to the integrated density of states of 
the Laplace operator on the product of a torus with Euclidean space. The 
asymptotics we obtain are then used to compute the density of states of the 
magnetic Schrodinger operator. 


1. Introduction and Main results 

The first problem we are considering in this paper has several equivalent formu¬ 
lations. 

1.1. Number theoretic formulation. For p > 0 and k £ M d , let B(p\ k) be the 
ball of radius p centered at k. Let S'(p;k) be the number of integer points inside 
the disk B(p, k) C R 2 . The classical Gauss Circle Problem consists in estimating 
the remainder term 

(1.1) R(p;0)=S(p;0)-7rp 2 

Hardy and (Edmund) Landau have found lower bounds for this problem, while the 
current best upper bound is given by Huxley in |_4j- This problem has also been 
studied for balls of dimension higher than two, see e.g. [5], and it is well-known 
that averaging over the radius of the ball improves regularity of the remainder. 

In this paper, we consider a variation on this problem: we estimate the measure 
of the intersection of affine planes sitting on integer coordinates with balls of large 
radius in R d . More precisely, put 

(1.2) A k := Z k x R d ~ k C R d 

and let B d (p 1 k) be a ball in R d of radius p centred at k := (ki,k 2 ) £ R fc x R z , 
where k + l = d. Denote by S(p; ki; d, k) the ^-dimensional volume of the set 
B d {p, ki) nA fc . A simple observation shows that we have 

(1.3) S{p- 1 k 1 -d,k)=uji ^2 {P 1 - |7 ~ ki| 2 ) i/2 , 

7 ez fe 

| 7 —ki | <p 

where ujd is the volume of the unit ball in R d . One can see that the integral of 
R{p, k) over k 2 £ T l = R l /Z l , is the same as the remainder term 

(1.4) R := S(p-, kp, d, k) — uj d p d , 
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obtained from Equation (1.31. Our aim is to compute an estimate of R for large 
values of p. Before discussing the results, we would like to describe different for¬ 
mulations of this problem. 


1.2. First spectral theoretic formulation. Let 

(1.5) H = -A + V 

be a Schrodinger operator acting in with a smooth real-valued periodic potential 
V ; for simplicity we assume that the lattice of periods T = (27rZ) d , with dual lattice 
rt = Z d Denote the integrated density of states (IDS) of H by N( A) := 

It can be defined by the formula 

(1.6) N(\;H):= lim N(X ’? l) , 

L—foo L a 

where Hl is the restriction of H to the cube [0 ,L] d with appropriate self-adjoint 
boundary conditions and N( A, Hl) is the counting functions of the (discrete) eigen¬ 
values of H^. Note that this parameter A is related to the parameter p of the 
previous section by p = \/\. While this formulation of the IDS is important for 
Theorem |1.5[ for periodic V we use an useful equivalent definition. 

Following m, we express H as a direct integral 


r(B 

(1.7) H= H( k)dk, 

JT d 

Then, one can express in terms of the counting functions of the fibre 

operators H{ k): 


( 1 . 8 ) 


tf(A):= 7 / JV(A;JT(k))dk, 

J T d 


where IV(A, H( k)) is the eigenvalue counting function of H( k). Remarkably, despite 
the fact that the asymptotic behaviour of N(A,H(k)) for fixed k and A —> oo is 
very irregular (so that even the precise size of the remainder 

(1.9) R( A; k) := N( A, H( k)) - C d X d/2 


is unknown), integration over all quasimomenta k € T d := R. d /Z d makes things 
extremely regular, so that there exists a complete asymptotic expansion of N( A) in 
powers of A as A —> oo, M- Here, we have denoted 


( 1 . 10 ) 


C d 


u d 

(27 l) d 


and oj d 


n d / 2 

T{l + d/2) 


is the volume of the unit ball in R d . The question we want to study is what would 
happen if, instead of integrating against all quasimomenta, we integrate over a 
subset of them, say over an affine plane. We write k = (ki,k 2 ), where ki £ T fc , 
k 2 € T* and define the partial density of states (PDS) as 


(1.11) N p ( A;ki) = AT p (A;ki;d, k) := . \ , f N(X, H(k 1 ,k 2 )) c\k 2 . 

(Zn) Jjl 

Our aim is to investigate the asymptotic behaviour of the PDS as A — > oo. 
Obviously, the regularity at infinity will be improving as l increases and so the 
larger l is, the more asymptotic terms we are likely to obtain. This asymptotic 
problem can be treated in two steps: 








DENSITY OF STATES 


3 


Step 1. Obtain the asymptotic behaviour of the PDS for unperturbed operator 
H° := —A. More precisely, we want to obtain as good an estimate on 

(1.12) i?°(A;ki -d,k) := IV°(A; kn; d, k) - C d X d/2 

as possible (of course, superscript 0 refers to the fact that we are dealing with the 
case V = 0). A simple calculation shows that if k = 0, then i?°(A; ki; d, 0) = 0, 
so this step is trivial when dealing with the IDS. In the case of k > 0 this step 
becomes quite non-trivial and interesting. Once we have performed this step, we 
can move to 

Step 2. Compute (or estimate) the difference 

(1.13) N P (X; ki; d, k) - N°(A; ki; d, k) 

and try to obtain as many asymptotic terms of it as possible. It follows from a 
simple computation that 

(1.14) A°(A; k i; d, k ) = {2ir)~ d S{V\; k i; d, k), 

hence the main aim of this paper deals with the first step of this programme; we 
intend to perform the second step in a separate publication. 


1.3. Second spectral theoretic formulation. Consider the operator H = — A + 
V acting on T l x with a smooth potential V : T* x R fc — > R. We assume that, as 
a function on V is periodic with the lattice of periods (2-7rZ) fc . Then, we have, 
from the definition of both the IDS and the PDS that 

(1.15) N( A; H) = [ N( A; H{ k 2 )) dk 2 = (2n) l N p (\- 0; d, k), 

\ 17 rj J t ! 


that is to say that the integrated density of states equals the partial density of states 
up to a constant. If we consider a more general (but also less natural) operator H ^, 
the domain of which consists of functions on T* x which become periodic after 
multiplication by e iklXl , then the IDS of H equals, again up to the same constant, 
7Vp(A;k i;d,k). We would also like to mention that expression (1.31 appears in the 
study of integer points in anisotropically expanding domains. This has applications 
in the study of the asymptotic behaviour of the eigenvalue of the Laplace operator 
on the torus in the adiabatic limit, and was developed in |5I. 


1.4. Main results. Our first main result is as follows: 

Theorem 1.1. The error term R(p', ki; d, k) satisfies the asymptotic estimates 

{ O (p^- 1 )/ 2 ) ifk<(d+ l)/2, 

(1.16) R(p-,ki;d,k) = < O (p (d-1)/2 logp) ifk=(d+ l)/2, 

I O (^p d ~i-d+2k j ifk>(d+ l)/2 

uniformly in kj. 

Remark 1.2. Recall that R(p; ki; d, 0) = 0 for all values of p, k 1; d. 

We do not pretend that all of these estimates are optimal, but some of them are, 
as can be seen from the following result: 
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Theorem 1.3. For k > 1 and p sufficiently large, there exists a positive constant 
C dyk and € T fc such that 


(1.17) 


R(p\ ki; d, k) > 


d-± 

F d ,kP 2 

n — 
Cd,kP 2 


if d = 1 mod 4 
else, 


where e > 0 is arbitrary. When d ^ 1 mod 4, £/ie lower bound i?(p;ki;d, fc) > 
Cd,kP^~ holds for k = 1. 


In particular, this theorem means that for 1 < k < and d ^ 1 mod 4, we 
cannot get improvements on the upper bounds found in Theorem |1.1| It also means 
that for d = 1 mod 4, /c ^ 1, we cannot get improvements in the exponent. 


Remark 1.4. It seems interesting that, after we have integrated N(A',H(\t)) (d — 
l)/2 times, additional integrations do not improve the remainder estimate, until we 
perform the last (d-th) integration, which makes the remainder equal zero. 

Open problem. The results in [2] imply that for k = d, our upper bound is not 
optimal, but as d -A oo, our upper bound converges to the optimal one, in the sense 
that d— (d — l _ 2 d+ ^ k ) —> 2. Hence we may ask what is the optimal upper bound for 
k > ^±1. 

1.5. Operators with constant magnetic field. Another type of problems we 
consider in this paper is the asymptotic behaviour of the density of states of the 
(Lev) Landau Hamiltonian (Schrodinger operator with constant magnetic field). 

Let Dj = —i-ggr- Then we define the Landau Hamiltonian H d as the operator 
acting in whose action is given by: 

Hd = (D\ + X 2) 2 + D\ + • • • + D d . 

Of course, only operators H 2 and H 3 make real physical sense, but for the sake of 
completeness we will deal with all dimensions. 

Let kl d {p) for d > 2 be the parabolic domain in given by 

(1.18) ^ d {p) '■= {(^Oi x) £ R d : 0 < x 0 < p — |x| 2 } . 

Defining P(p; d, k) analogously to S(p ; 0; d, k), that is, 

(1-19) P( P ;d,k) = Vol l (n d (p)nA k ), 

one can see that 

LpJ 

(1.20) P(p; d,k) = Y, S({p - j) 1/2 ; 0; d - 1, k - 1). 

j =0 

The IDS N( A; H d ) is related to P(p; d , k) by the following proposition. 


Proposition 1.5. Let H d be the d-dimensional Landau Hamiltonian. Then, its 
integrated density of states is given by 

( 1 . 21 ) N(\-,H d ) = 2^n 1 - d p(^±-,d-l,l) 


for p > 1 , and 0 otherwise. 
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We get an asymptotic expression for P{p\d,k), via the next theorem. Defining 
E 0 (p) := Eo(p,d) = ^ lP {d+1)/2 + ±P id ~ 1)/2 and 


E n {p) := E n (p, d) = E 0 + 

k =1 


s 2 k r(^±i) 

(2fc!) r(^±^SE) P 


we obtain the following theorem. 


Theorem 1.6. As p —> oo, P(p;d,k) admits the asymptotic expansions: 

(1.22) P(p; d, 1) = Ud-iE^d+i^ (p) + 0(1), 


(1.23) = 

If k > , we have 

(1.24) P{p- d, k) = E, |*» , (p) + OipH^-mm)). 

Finally, if k < 

d+4 r 

(1.25) P(p;d,k) = E^^p) + 0 {p~ (log p) 5 ), 

where <5 = 1 if k = and 0 otherwise. 

Replacing the result in Proposition |1.5| with the asymptotics in Theorem |1.6| we 
immediately deduce the following corollary. 

Corollary 1.7. The integrated density of states of the Landau Hamiltonian on M 3 
admits the asymptotic expansion 

N(X-H 3 ) = ^X 3 / 2 + 0(l) 

for large enough X. 


The rest of the paper is organised as follows: in Section 2 we formulate several 
results which will be used in the proof of the main theorems, but we will postpone 
their proofs until Section 6. In Section 3 we prove the upper bounds in the Laplace 
case, and in Section 4 we obtain lower bounds. Finally, in Section 5 we deal with 
the magnetic case. 
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2. Auxiliary results 

In order to prove Theorem 0 it will be useful to give an alternate expression 
for S(p; ki; d , k). Let us define the function x : R fc —>■ M as 


( 2 . 1 ) 


x(x) = 


(1 — \x\ 2 ) 1 / 2 if |cc| < 1, 


0 otherwise. 

We can then observe that 

(2.2) S(p- 0; d, k) = uip 1 Y x(n/p). 

n£Z k 

We would like to use Poisson’s summation formula 

( 2 -3) Y /( n ) = Y /( m ) 

n&L k meZ fc 

with / = x- This will allow us to get upper bounds for all ki £ T fc , from the 
relation 


(2.4) 


T(/(a; - ki)) = e 


_ 27riki-£ 


H*M), 


where T is the Fourier transform operator. For the rest of this section, we therefore 
consider ki = 0, and it will be seen in the proof of Lemma |2.2|t hat this assumption 
is made without loss of generality. In order for Equation (2.31 to hold we need to 


smooth out x- To do so, we will consider its convolution with Friederichs’ modifier 
Hence, setting Xe = 4' e * x we get that 


(2.5) 


Xe(0 = -Mom 


Theorem 1 1.1 1 follows from two lemmas. The first one finds asymptotic upper and 
lower bounds for S: 

Lemma 2.1. Let xt and x7 be defined on M. k by 

(2-6) xt( x ) = ^^ e y Xe((l T e)x). 

Then, we have that 

(2-7) X7 ( x ) < X(x) < Xt( x ) 

for all x £ R fc . Immediately, if we define 


( 2 . 8 ) 


we get that 


St(p) =UJ iY xt(n/p), 

ne Z fc 


(2.9) S~(p) < S(p) < 3+(p). 

Since Xe are smooth functions, we can use Poisson’s summation formula to 
compute the asymptotic expansion of Sf 1 . The second lemma therefore gives the 
asymptotic expansion of x(£). 
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Lemma 2.2. The Fourier transform of x satisfies 

< 2 ’“) ««) = cos ( 2 ' KI * + °(l«r l ' ,+3,/! ) 

for some C > 0 as |£| —> oo. Furthermore, its derivative satisfies 

* 2,11 * 4l« {) = si ” ( 2,r|?l “ f) + °(ir 

In particular, the asymptotic behaviour of both x(f) and its derivative does not 
depend on the co-dimension k. 

We will postpone the proof of these lemmas until Section [6] 


3. Proof of Theorem 11.1 


In this section, we prove Theorem 1 1.1 1 using both Lemmas |2.2| and |2.1| We have 
that 

(3-1) S~(p) < S(p) < S+(p). 

Let us therefore find asymptotic expansions on Sf 1 . We shall split those computa¬ 
tions in two cases : whether k > (d + l)/2 or k < (d + l)/2 

3.1. Case 1. Here, we assume that k > (d + l)/2. Let us find asymptotic expan¬ 
sions on Sf 1 . Since Xe is a smooth compactly supported function of x , we may use 
Poisson’s summation formula (|2.3|) to obtain 


(3-2) Sf 1 = to d p l Y xt( n /p) = Lo d p d Y xt(pm). 

n£Z fc m£Z fc 

Since we have that 


(3-3) 


1 


Xt(mp) = {1 ~ e)d *(zmp)x({f- e ) 


we get, assuming e 1/p, that 
(3-4) 

Sf =u d Y ( 1 + 0(e))p d T(emp)x(mp) + 0 | Y ernp d+1 ^{emp) \x'(mp)\ ] , 

me 1 k \m€Z k 

which directly implies 


( 


Sf = Lo dP d + 0{ep d ) + O 


(3-5) 


\ 


Y P d ^( em P) \x(mp)\ 


me Z fc 

\|m|^0 


/ 


+ 0 | Y emp d+1 'I>(emp)\x , {mp)\ ) • 

\mez fc ) 

Observe that T(£ ) = Od^l 9 ) for any q whenever |£| > 1 and bounded for |£| < 1. 

that %(£) = 0(|Cr (d+1)/2 )- Hence, choosing q = d_2fe_1 


Recall from Lemma 


2.2 


2 > 
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the third summand in (3.5) can be split into two terms, becoming 
(3.6) O 


/ 


1 i V 


\ 

n (d~ l )/2 

v 

1 


P 

z^ 

1777,1 (d+l )/2 

( ep )( 2 k+ 1 -d)/ 2 | m |fc + 1 

J 

V 

_l<|m|<l/ep 

|m|>l/ep 



The first sum can be estimated by 

1 f 


E 


!/ e P r k -1 


(3.7) 


l<|m|<l/ep 


l|(d+l )/2 


r (d+ 1)/2 


dr 


|o(m + 

\o(log ep) 


if k > (d + l)/ 2 , 
if k = [d + l)/ 2 . 


The second sum can be estimated by 

1 


(3.8) 


E 

mez fc 

|m|>l/ep 


( £j 0 )( 2 fc-d+l)/ 2 | TO |fc + 1 


„k-l 


/l/ 6p ( e/ 9 )( 2 fc-d+l )/ 2 r fc+l 


/, , d+l- 2 /c\ 

dr = O ( (ep) 2 J . 


As for the last summand, it is easy to see with the same computations and using 
4*(£) = O ^|£| 2 ^ that the extra power of ep exactly compensates the extra 

power of m, and we have that the asymptotic behavior in ep is the same for all 
for summands whenever k > (d + l)/2. Furthermore, when equality holds, the 
polynomial component is the same. Therefore, we have to choose e = p~ 3 such that 

(3.9) " 

This is achieved exactly when 

(3.10) j = 


e P d = pl d -W(e P y 


2 k 


1 — d + 2k 

This gives us the announced asymptotic estimates when k > (d + l)/2, that is 

u>dp d + 0 (p d ~ ) if k > (d+ l)/ 2 , 

UdP d + 0(p^ logp) ii k = (d+ 1 ). 


(3.11) 


S(p) = 


3.2. Case 2. We now assume that k < (d + l)/2. In this case, we have that the 
sum converges with T = 0(1). Hence, the asymptotic expansion for Sf 1 simplifies 
to 


(3.12) 

Sf = L0 d p d +O{ep d )+O 


( 



( 

- \ 

n (d~ l )/2 V 

l 

i n 

jd-D/2 sr 

ep^(emp) 

p 

| m |(d + 1 )/2 


p zL 

i TO |(d-i )/ 2 

raEZ fc 

\ \m\^0 

) 


rnEZ fc 
\ ra^O 

/ 


The third sum converges and the last one as well if k < In that case, choosing 
e = p-O+i )/ 2 satisfies Theorem and choosing e smaller does not improve the 
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estimate. If k = | or k = ^7 using 4*(£) = O (|£| _1 ) for m > (ep ) _1 yields the 
same result, finishing the proof. 

Note that Equation (2.4) ensures that these estimates hold for all ki £ T k . 

4 . Lower bounds 

Let us first follow the argument given in [I] for d = k = 2. The beginning of the 
argument is the same, which we add for completeness. Since R(p\ ki) is periodic in 
ki with respect to T, we can compute its Fourier coefficients, obtaining 
(4.1) 


[ R(p-k^e- 2 ^^ dk, = f [ 
J T fc JT k l 

[ P l x(—)e- 2 ^^ dk, 

jR k V P J 


-W d p d + p l ^2x 

7GT 


7 - kj 


- 1 e _27rikl ' 7 dk. 


= P 


TTJT 1W9 cos ( 2 ttp| 7 | - ^ d+1 '> 7r \ + 0 (| p7 | (<*+ 3 )/ 2 ) 


,(p|7l) (d+1)/2 


from Lemma [ 2 . 2 | Additionally, we have that 


(4.2) 


[ R(p; ki) dki = 0 . 

JT k 


Hence, for all 7 £ T \ {0}, we have that 
(4.3) 


J T k 
> max 


I R(p; ki)|dki 

( [ R(p> k 1 )e _27rikl ' 7 dk, , [ R{p\ k 1 )e _4,rikl ' 7 dki 
\ Jj k J T fc 

cos (2np\j\ - 1 ) 7r > j ) _T r cos(dnp\^\- 


P 2 

> C max 


(d + 1)7 


7 2 

for C, c positive constants whose value can change throughout. Whenever d ^ 1 
mod 4, we have that 


(4.4) 0 < inf max 


cos x — 


(d + 1)7 


cos 2 x — 


(d + 1)7 


hence in that case, fixing 7 £ T, we conclude that there exists r* such that for all 
r > r* 


(4.5) 


/ T k 


\R(p', ki) | dki > Cp 2 


We conclude that whenever d ^ 1 mod 4, 


(4.6) 


sup I?(p; ki) > Cp 2 

kiGT fc 


The remaining case, that is when d = 1 mod 4 is more subtle. We will use 
results found in [7][Theorem 3.1, Lemma 3.3]. Indeed, from Equation (4.31, we 
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have 

(4.7) 


/ \R(p; kj.) | dky > 

Jjk J 2 


COS 


(27rp|7| 


7" 


From Lemma 3.3 in [7], we know that, if k > 2, for all e > 0, there exists po > 0 and 
a £ ( 0 , 1 / 2 ) such that for all p > po there exists 7 £ F such that I 7 I < (2np) e and 
the distance from 2pj to an integer is greater than a. Choosing such a 7 bounds 
cos(27rp|7| — tt/2 ) away from 0 , and we get that 


(4.8) 


/ |/2(p;k 1 )|dk 1 >C'p i i 

Jf k 


d — 2 d+1 


Since e > 0 is arbitrary, we get the desired result. 


5. An APPLICATION TO THE LANDAU HAMILTONIAN 

5.1. The Landau Hamiltonian. Decomposing H d = H 2 CD Dd- 2 , we can first 
study the problem 

H 2 u = Xu. 


Consider the definition (1.6 1 for N(X;Hd), with periodic boundary conditions for 
X\ and Dirichlet boundary conditions for x = (x 2 , • • •, x d ). 

For H 2l we can write the solutions as u(x i,x 2 ) = e~^~ Xl f(x 2 ), which reduces 
the problem to solving the eigenvalue problem 

((Ci + x 2 ) 2 + D%) f(x 2 ) = A/(x 2 ). 

This is a shifted quantum harmonic oscillator. We have that cr(H 2 ) = {2j + 1 : j £ 
N}, each with infinite multiplicity. It is a standard computation, see e.g. [5], that 


(5.1) 


N(X-,H 2 )=- 


A- 1 


for A > 1, and 0 otherwise. Extending the methods of [ 6 ] to higher dimensions, it 
is again a simple computation to show that for A > 1, 


(5.2) 


N(X-,H d ) = 


Ud -2 


L^J 


( 2 tt)< 


( A “ 2n ~ !) 


(d- 2)/2 


n —0 


Thus, from the definition of P(p', d, k), we have indeed that 


(5.3) 


N(X-H d ) = 2^7T 1 


A - 1 


~\d — 1,1 


5.2. Computations for general paraboloids. In this section we prove Theorem 
|1.6| Consider the expression 

LpJ 

(5.4) P(p; d, k ) = £ S((p - j) 1 / 2 ;0; d - 1, k - 1). 

i=i 

By Theorem 0 we have 

LpJ LpJ 

Yl S ((P-ti 1/2 ’ 0 ’d- hk-l) = Y^ (vd-i(p-j) (d ~ 1]/2 + 0(X(p))J , 


j=0 


j=0 
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where 


(5.5) 


X(p) = 


'pl( d 1 2 k—d ) if k > (d + 2)/2, 
p( d ~2)/4 i 0 g p Hk = (d + 2)/2, 
j{d-2)/4 if 1 <k< {d + 2)/2, 


if k = 1. 

Comparing with the integral, we get that for all X as defined above, 


(5.6) 


LpJ 

Y J X(p) = 0(pX(p)). 

3=0 


For any d, we can use the Euler-Maclaurin formula : 

6 f(a) + f(b) 


(5.7) 


E /( n ) = [ dx + 

^ J a 


E 

k=l 


& 2 k 

W 


j2/c—1 . 


dx 


j2fc—1 . 


x—b 


dx 


O 


d 2p f 


dx 2p 


d t 


for any integer p > 1, where Bk is the fcth Bernoulli number. Note that for integer 
a, 

(5.8) = £/--'>/*. 

3=0 i =0 

Hence, by the Euler-Maclaurin formula, we get that 


(5.9) 


E(«^) (d-1)/2 

3=0 

= r^- i )/ 2 dt+ 

Jo 


(d— 1)/2 


y- JkjEE a ^ +0 ( a -i/2) 

/ ^ /o7„!\ TD/d+3-4fc\ a -r /• 




' (2fc!) r( 


Obviously, when d is odd, this last sum is actually finite and the error term 0. 

When p is not an integer, we can write p = a + r, where r is the fractional part. 
In that case, using the Euler-Maclaurin formula again, we get 


B «+-- i ) w - 1)/2 = Bi + r >"‘ _1)/2 

3=0 3 =0 

= J\t + t )^- 1 )/ 2 dt + \ ( r ^- 1 )/ 2 + p ^- 1 )/ 2 ) 

B 2k IW) 


+ E 


i ( 2 fc!) r( d+3 ~ 4fc ) 


/ d+1 


2 — 7- 2 


+ 0{t) 


d 

+ E 


- ^pC +1 )/ 2 - r (rf+ 1 )/2^ + £ ^ r (d-l)/2 +p (d-l)/2^ 

r(^±i) 


fc<- 


t ( 2 fc!) r( d+3 ~ 4fc ) 


/ d+i^ 


2 — 7- 2 


+ 0(r). 
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Let us observe that 

^p ( rf+ 1)/2 + E“= 0 (P-j) (d ~ 1)/2 


lim 

p—>oo 


= lim 

p—too 

= 1. 


l p ( d ~l )/2 

- - —/y-(rf+l)/2 + T (d- 1)/2 + p (d-l)/2 + 0(p(d- 3)/2) 


P' 


4-1 


This is because r = 0(1). Similarly, if we define P 0 = OL-p( <i + 1 )/ 2 + ^ p <d L/ 2 and 


d+i 1 

F = F | y r (^) ^|r 4 » 

n 0 ^ ( 2 fc!) r( d+3 ~ 4fc ) P ’ 


we get that 


lim 

p—>• oo 


-^ + E“= 0 (p-i ) (<i - 1)/2 


p(d-l)/2-2n-l 


^(n+i) 1^ 

(2(n + 1))! 2 j2n+1 


whenever (d —1)/2 — 2 n — 1 > 0, after which point the contribution of the fractional 
remainder r gets more important than the denominator. Hence, we obtain the 
asymptotic expansion 


LpJ 




(d— 1)/2 _ 


1)/2 . 


1 


(5.10) 


i=o 


d + 1 


;P 


(d-l)/2 


B'2k r( 2 ) . 

(2k\)T{ d + 3 - 4k ) P 2 


When k = 1, we already have that X(p) = 0. Therefore, we have that 
P(p, d, 1) 2 mxi \ /o 1 


Ud -1 


_ ft ^ 1 )/ 2 , L n ( d - 1)/ 2 

d + r ^2^ 


+ E 


p. 


2fc 


T(^) 


l<fe<- 


_ 3 (2fc!) p(d+T=4fc) 


d+1—4fc , s 

P 2 +0(r), 


from which we recover a (quite sharp) asymptotic integrated density of states for 
the magnetic Hamiltonian H^+i- 

Let us combine equations ( |5.5[ ) and (5.10). When k = d, we get that the error 
term from X is greater than E+ and as such, 


P(p; d, d) = 2 j d 1 + O 
d+1 


( d 2 -d+ 2 \ 

P 2d J • 


When fc > EE we get that 


P(p; d, fc) 


2 d+1 1 d — 1 

P 2 + ^P 2 


d + 1 

+ E 


P 


21 


r(E 1 ) 


i <j<£? 


_ t (2j)! r(^^±E iI ) 


9 I d+1 — 4j 

347TP a 


+ 0(pE d 1 2fc-d)). 


<+2-1 
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Finally, when k < ^± 2 , we get that 

P(p ; d, k) 2wd_i a±i 1 <1=1 
-— = j , i P 2 + oP 2 

Ud—i d -hi 2 


E 

i <o<H 


B 


2 j 


( 2 j)! r(d±| = 4i 


d~l~ 1 \ 

- d+1—4j , d+4 , , 

YZ7T P 2 +0{p^r (log p) S ), 


r(T) 


where <5 = 1 if k = 4±2 and 0 otherwise. 


6. Proofs of auxiliary results 


6.1. Smoothing of the cut-off function. Let us define a smooth, even bump 
function ip in C^°(K), supported in [—1,1], such that the integral 

r°° i 

(6.1) / ip{r)r k ~ x Ar = -—, 

Jo Vk -1 

where 14_i is the area of the unit sphere in R fc . 

Using this function, we can define the radial bump function il f e on K fc , of total 
mass 1 to be given by 

(6.2) * e (x) = 

e K e 

Let iff := i and Xe( x ) = \El(x) * xi x )- Its Fourier transform is given by 


(6.3) Xe(0 = *&m)- 

Let xt and x7 be defined on K fc by 


(6-4) xf{x) = , * v Xe((l -F e)x). 

(1 =F e)' 

We can now proceed with the proof of Lemma |2.1| 


Proof. To show that Xe ( x ) — xi x ) — xt > the idea is to obtain xt ( x ) by averaging 
x(x) on a ball of radius 0 < e < x about each x. To do so, first notice that 


(6.5) 


Xe(z) < sup (x(x-t)) V e (x)dx 

|t|<£ J R fc 


if \x\ < e, 

(1 — (\x\ — e) 2 )s if e < |rr | < 1 + e. 


If we show that 

( 6 . 6 ) 


Xe{x) < (l + e) l x 



we get the desired lower bound. Indeed, taking y = in the preceding equation 
yields 

(6-7) x(y) > n * w Xe((i + Oy) = x7(v)- 

(1 + e ) i 

Therefore, it only remains to show that ( |6.6| ) holds for all x £ R fc . First note that 
if 1 2 1 > 1 + e, both sides are 0. We shall split the remaining cases in |x| < e and 
e < |a;| < 1 + e. 
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Restricting ourselves to the first case, if \x\ = e, we get that 


( 6 . 8 ) 



= (l + 0‘ 



= (l + 2e)5 
> 1 


1 

2 


> Xe ( x ) 


Since x(yiy) is a decreasing function of |x|, we conclude that ( 6.61 holds for 0 < 

M < e- 

In the case where e < |cc| < 1 + e, we need to show that 


(6-9) (1 - (|ar| - e) 2 )^ < (1 + e? (l - “ . 

It is equivalent to show that 1 — (|x| — e) 2 < (1 + e) 2 — |x| 2 . This is the case if 

1 — \x\ 2 + 2\x\e — e 2 < 1 + 2e + e 2 — \x\ 2 
(6.10) <^>2|x|e< 2e(l + e) 

<t=> |x| < 1 + e. 


Since the last line is true by hypothesis, we can conclude that the left-hand side 
inequality of (2.71 is true. 

In order to get an upper bound on \(x ), we proceed in a similar fashion, averaging 
Xe{x) on a ball of radius e around x, which yields 


( 6 . 11 ) 


Xe(x) > inf x(x - t) 
\t\<e 


> 


(1 -{\x\+e) 2 Y 
0 


if |x| < 1 — e, 
otherwise. 


As we did before, it suffices to show that 

(6-12) Xe(x) > (1 — e)*x (j~) • 


Notice that the left hand side of that equation is 0 whenever |x| > 1 — e. Like 
before, we see that 


(6.13) 


(l-(|x|+e) 2 )*>(l 


eY 



1 

2 


is equivalent to \x\ < 1 — e. This concludes the proof. 


□ 


6.2. Fourier transform of y. 

Proof. Let us compute y(£). We will split the cases k = 1, k = 2, and k > 2. If 
k = 1, then 
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X(0 = J 1 (1 - a; 2 )( d -i)/ 2 e- l2TO « dx 


(6.14) 


C 

W 


2 J d/ 2(2tt|£|) 


C 


cos ^27 t|£| — ) + 0(|C] (d+3)/2 ) ; 


|£|(d+l)/2 

using 0[Eq.3 .387 and 8.451], which is the desired result. 
We also obtain that, following [3][Eq. 3.621] 


(6.15) 


X(0)=2“B(i±l 1±!). 


2 2 

Using identities of the Gamma function, we get that 


(6.16) 


^(^±l^±i) = 


T d /2 


r(f+ i) 


= Wd, 


which is the desired value. 

If k = 2, then the Fourier transform is given by 

(6.17) X0£)= [ X(x)e- i2nx< dx. 

J R 2 

Working in polar coordinates, we get that 


(6.18) 


n Z7T 

r(l 

= [ r( 1 - r 2 


_ r 2^(d-2)/2 e -i27rr|^| cos 0 

)( d - 2 )/ 2 j 0 (27r|£| r ) dr 


d0 dr 


|£|d/2 ^d/2(2 7r l^l) 

= |^|(d+l)/2 C0S (Wl - + 


which is the desired result. 0[Eq. 8 .411, 6.567 and 8.451] were used respectively 
for an integral formula for the Bessel function, its integral, and its asymptotic 
expansion. 

We also obtain that 

27r 

(6.19) X(°) = T - 

Using identities of the Gamma function, we get that 


( 6 . 20 ) 


27T 7T d / 2 

Ul a ~ r(| + i) 


which is the desired value. Finally, if k > 2, then, working in spherical coordinates, 
we get that the Fourier transform of x is, for some constant C, 
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x(0 


cf f r fe - 1 (l-r 2 )' /2 sin fc - 2 6>e- i27rr| « |cose d6»dr 

JO JO 


= |^|(fc- 2)/2 J Q rk/2{l ~ r2 ) l/2j (k- 2)/ 2 (27r|e|r) dr 

21 ) C ! 

= |£|(fc-2)/2 |£|(Z+2)/2 Jd/2(27 r |^'|) 

= |g|(d+l)/2 C0S ( 27r l£l - + °d^ (d+3)/2 )- 

using j33 [Eq. 8.411] in the first line, which is the desired result. 
Additionnally, we have that 


( 6 . 22 ) 


X(0) = Vol(S* _1 ) [ r fe-1 (l — j- 2 )( d_fe )/ 2 dr 

Jo 

_ 7T fc / 2 B(|, ^f±^) 


r(f) 


Using identities of the Gamma function, we get that 


(6.23) 


X(0 )ud-k = u d 


which is once again the desired value. 

One can note that in each of those cases, we ignored the trigonometric term to 
get an upper bound, considering it to be 1. Hence, since translation by ki is simply 


multiplication by a complex exponential in Equation (2.3), it can be ignored in just 
the same fashion. 

Finally, we get the result for the derivative using the identity J' u = J v -\ — Jv+i) 


and basic trigonometric identities. This completes the proof of Lemma|2.2 


□ 
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